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In order to describe a nonuniform equilibrium mixture with an interface between two coexisting 
phases it is necessary to consider contributions to the Helmholtz energy which depend on the gradi- 
ents of for instance the density. Van der Waals [l], @] was the first to introduce such a term, which is 
very important in the interfacial region for a one-component system. Cahn & Hilliard [j§] extended 
£ — ■ this analysis to a binary mixture by introducing gradient terms of the mol fraction. We give an 

systematic extension of the gradient theory to three-dimensional multi-component systems. 
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INTRODUCTION. 



In order to describe the equilibrium properties of an interface between two coexisting phases, using a continuous 
model, it is necessary to consider contributions to the Helmholtz energy which depend on the gradients of for instance 
the density Van der Waals [l|,[3] was in 1893 the first to introduce such a term for a one-component system. In 
1958 Cahn & Hilliard Q extended the analysis of van der Waals and introduced gradient terms of the mol fraction 
<L|__i in binary mixtures. As the Helmholtz energy density given by van der Waals is no longer a function of the local 
density or local densities alone, there is no local equilibrium in the traditional sense in the interfacial region. The 
continuous description is in other words " not autonomous" . We refer to the monograph by Rowlinson and Widom [3| 
for a thorough discussion of the van der Waals model in general and of this point in particular. A lot of work on the 
equilibrium gradient model was done by Cornelisse [5J. We refer to his thesis for the relevant references. 

The gradient model is often used for a system, in which properties vary only in one direction. We do not restrict 
ourself in this manner and do the analysis for the three-dimensional system. 

We use the standard thermodynamic variables which obey the standard thermodynamic relations for homogeneous 
mixtures. In Sec. [II] we consider the inhomogeneous mixture. We postulate the dependence of the specific Helmholtz 
i energy on the thermodynamic variables and their gradients, using the fact, that in equilibrium the temperature is 

constant. Using that in equilibrium the total Helmholtz energy has a minimum and that the amount of the various 
components is fixed, expressions are derived for the chemical potentials of all the components, using Lagrange's 
method. This is done for a different choice of variables in Sec. [Ill j and Sec. [IVj . Extending a method developed 
by Yang et. al. @ an expression is also found for the pressure tensor. We derive the different forms of the Gibbs 
relations. Their importance is crucial for non-equilibrium description which should be based on the equilibrium 
analysis. Explicit expressions are also given for the internal energy, enthalpy and Gibbs energy densities. In Sec. |VI| 
we give a discussion and conclusion. 

In following papers we aim to generalize the analysis to non-equilibrium systems, using the results of this paper. 
This will extend the work of Bedeaux et. al. 0, [H, Q for one-component systems, in which the properties varied only 



in one direction. 



II. THE GRADIENT MODEL. 



In order to describe inhomogeneous systems in equilibrium, one could assume, that this can be done by the usual 
thermodynamic variables, which depend on the spatial coordinates. All standard thermodynamic relations are then 
assumed to remain valid. 

As van der Waals has shown for one-component system however, this is not enough to describe the surface. It 
is necessary to assume that thermodynamic potentials, particularly the Helmholtz energy density, also depends on 
spatial derivatives of the density. Cahn & Hilliard have shown for a binary mixture [3j, that the Helmholtz energy 
should depend on the gradients of the mole fraction of one of the components. For a multi-component non-polarizable 
mixture we shall use as general form of the Helmholtz energy 

1 - 

/(r) = / (T; zr(r), . . . , z n (r)) + - £ ««(zi(r), . . . , z n {v)) V*(r) • Vz,-(r) (II.l) 

«ti=i 

where /o is the homogeneous Helmholtz energy and Zfc, k = l,n, where l,n indicates all integers from 1 to n, are 
generalized densities. These can be either total concentration of the mixture together with the n — 1 fractions of the 
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components, {c, £1, . . . , (see Sec. or n concentrations of the components, {ci, . . . , c n }, (see Sec. 

Concentrations may be taken either on the molar basis or on the mass basis. 

This form of the specific Helmholtz energy is obtained by the following arguments (cf. also 0, 0, E|). Dependence 
of the specific Helmholtz energy on the density gradients can be represented by Taylor series in these gradients. We 
describe here isotropic fluids, so any coefficient in this Taylor series can not depend on any direction and thus should 
be scalar. The zeroth term, taken when all the gradients are equal to zero, is the homogeneous Helmholtz energy 
/o(T; zi(r), . . . , z n (r)). In equilibrium, the total Helmholtz energy of the system has a minimum. Thus, the first 
order term, with the first order density gradients, V£j(r), is zero. The second order term is quadratic in the first 
order density gradients Vzj (r) • Wzj (r) , and is linear in the second order density gradients V 2 Zi(r). The latter one, 
however, contributes to the total Helmholtz energy the same way as the former one: 

J v cLyki V 2 Zi(r) = - J v dr Vac*- Vzj(r) + j v dv V-(k,- V^(r)) = 

n _ (II.2) 

= - Iv dr £ ir Vz -j ( r ) • Vz * ( r ) + Is dS K * n « ■ Vz * ( r ) 

i=i 3 

The first term on the right hand side of this equation can be combined with the quadratic term in the first order 
gradients. The second one can be chosen equal to zero by proper choice of the boundaries of integration. Thus, we 
end up with Eq. (jll.ip for the specific Helmholtz energy where the coefficients are the combinations of those from 
the quadratic in the first order density gradients term and the corresponding ones from the linear in the second order 
density gradients term. For ease of notation we will write Kij instead of Kj,-(Si(r), . . . , z n (r)), remembering their 
dependence on these variables. 

We note the ambiguity in the definition of the specific Helmholtz energy. Different expressions for /(r) give the 
same expression for the total Helmholtz energy F, due to the cancellation of the boundary contributions. This can 
be interpreted such that the measurable quantity is only the total Helmholtz energy, but not the specific one. To 
build the local description we need the local quantities, however. We will use Eq. for the specific Helmholtz 

energy, remembering that a divergence of a vector field, the normal component of which is zero on the boundary, can 
in principal be added. We will return to this point in the Appendix [X] . 

We may choose the matrix to be symmetric with respect to the component number (k^ = without loss of 
generality (since it appears only in symmetric combinations). We shall always take independent of the temperature. 

We note, that the gradient model is, as it has been used here, a general approach and is not only applied to surfaces. 
It has, for instance, been used in the description of critical behavior using renormalization group theory [Tol. Hlj|. In 
this paper we will focus on its use for the description of the surface. 

In the following paper, where we extend the analysis to non-equilibrium two-phase mixtures, we need all the 
thermodynamic variables. We derive all thermodynamic quantities and relations for the given choice of the independent 
variables. This is done in Sec. |III| for molar specific variables. We determine how the Helmholtz energy varies with 
a change of the variables and with a change of position and obtain so-called Gibbs relations. In Subsec. [Ill Cj we 
determine the physical meaning of the Lagrange multipliers and other quantities, for which expressions were derived. 
In Sec. |IVj we derive the results for volume specific variables following the same procedure as in Sec. [HI] . 

III. GRADIENT MODEL FOR THE MOLAR VARIABLES. 
A. The Lagrange method. 

We write the Helmholtz energy as 

/(r) = f (T, c, + K{c, & Vc, V0 (HI.l) 

where 

1 n— 1 1 n—1 

JC(c, f , Vc, V£) ee tl I Vc | 2 + - Vc- V& + i Yl — V &- V & ( IIL2 ) 
i_ ( c ^ c 

2—1 — 1 

Here and further we suppress r as an argument where this is not confusing. We also use £ as short notation instead 
of whole set {£i, . . . , and V£ as short notation instead of {V£i, . . . , V£„_i}. The molar density distributions 

are such that they minimize the total Helmholtz energy 



F= ( rfrc(r)/(r) 
Jv 



(III.3) 
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Assuming that no chemical reactions occur, the total number of moles of each component, vi = J v dr ^ (r) c(r) for 
i = l,n— 1, as well as the total number of moles, v — J v drc(r), are constant. The problem of minimizing the 
functional (|III.3|) . having n constraints can be done using the Lagrange method. Thus we minimize the integral 

„ n—l „ 

= / drc(r)[/(r)- Mn -y>i6(r)l =- / drp(r) (III.4) 
Jv i=1 J Jv 

where /i„ and ipi are scalar Lagrange multipliers. The concentration distributions which minimize the integral (|III.4|) 
must be solutions of the corresponding Euler-Lagrange equations. These relations give for the introduced Lagrange 
multipliers: 



n—l n—l 



(III.5) 

Q 1 n — 

^fe = (/o + /c) - - V- (n k Vc + ^ K ik V&) , & = l,n-l 



and an expression for p : 



p(r) = c 2 — (fo + IC) -cV-(«Vc + 5^«iVei) (IIL6) 



The ambiguity in /(r) discussed above does not affect the expressions for /i„, V'fc and p(r) (see Appendix [A]). We 
will see in Subsec. [Ill C| that the Lagrange multipliers fi n and ipk are the chemical potentials of the components and 
p(r) is a pressure. The exact meaning of this pressure as well as the meaning of the other quantities derived in this 
Subsection will be discussed in Subsec. |III C| . 

Solving Eq. (|III.5|) for c and £, one obtains the density profiles for the system. To do this one needs the values for 
the Lagrange multipliers fj, n and tf>k. 

Multiplying the first of the equations in Eq. (|III.5P with Vc(r) and the other (n — 1) ones with V£fc(r) and summing 
them all up we obtain the following expression 



(III.7) 



dx a 

where we use the summation convention over double Greek subscripts. The tensor 

<j a/3 (r) = p (r) 5 af3 + 7 Q/ 3 (r) (III.8) 
will be identified as the pressure tensor. Furthermore the tensor 

Wr) =j^^ + "£j?mm + + g K „ ?m?m ( „,. 9) 

OX a OX/3 V OX a OXp OX a OX/3 ' ~_ x OX a OX/3 

will be referred to as the tension tensor 2 . We note, that both a a f3(r) and 7^ (r) are symmetric tensors. 

From the definition (|III.4j) of the p(r), we can see that the Helmholtz energy given in (|III.lj) and quantities which 
are given by the Eq. 5[) and Eq. (|III.6|) are related in the following way 

n-l 

/(r) = Mn + J2 - P (r) »(r) (III. 10) 

i=l 

Including a possible divergence term in /(r) would modify this equation by adding this term also on the right hand 
side. 



1 This method is a generalization of the one, given by Yang et. al. [f| for a one-component system. 

2 The explicit expression for the pressure tensor in the square gradient model for a multi-component mixture was, to the best of our 
knowledge, not given before. 
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B. Gibbs relations. 

Consider the variation of the total Hclmholtz energy SF[T, c(r), £(r), Vc(r), V£(r)] with respect to the variation 
of the variables, which it depends on: 

SF[T, c, e, Vc, V£] = f dr {f(T, c, £, Vc, V£) 5c + c 5/(T, c, £, Vc, V£)} (HI.ll) 

where 

f /(T, c, fl Vc, VO^^iT+^fc+S^^ + ^iVc+j:^ *Vfc (IIL12) 

z=l 2=1 

is the total thermodynamic differential of the specific Hclmholtz energy with respect to the thermodynamic variables, 
which it depends on. Given Eq. Eq. pil.5|) and Eq. (|III.6|) . Eq. piI.12|) becomes 

*/(T, c, Vc, VC) = H <5T + 4 <5c + V yj k <% + (III. 13) 

where 

7i — l 71 — 1 n — 1 

50(c, £, Vc, VO = (k Vc + «* v &) ^ + H Vc + H K * v &) ( IIL14 ) 

i=l k=l i=l 

The total Hclmholtz energy variation becomes then 

n i~\ J> Tl 1 

SF[T, c, Vc, V£] = / dr {/(T, c, £, Vc, V£) 5c + c U± ST + J 5c + ]T V* } (m.15) 

c i=i 

since the boundary integral J s dS n s -S& disappears, because we have chosen boundaries of the system such that the 
density gradients are zero everywhere along the boundaries. Thus, we will interpret the expression in parenthesis as 
the total thermodynamic differential of the specific Helmholtz energy: 

r, r n — 1 

Sf(T, c, £, Vc, V£) = <*T + J Jc + £>fc <% (HI.16) 

2=1 

We note the ambiguity in the definition of the total thermodynamic differential of the specific Helmholtz energy. 
Different expressions Eq. (|III.13[) and Eq. (|III.16|) for Sf give the same expression Eq. 15|) for SF, due to the 
cancellation of the boundary contributions. This can be interpreted such that the measurable quantity is only the 
total thermodynamic differential of the total Helmholtz energy, but not the total thermodynamic differential of the 
specific one. This ambiguity is similar to the ambiguity in definition of the specific Helmholtz energy. We will use 
Eq. (|III.16[) . remembering this ambiguity. 

We write Eq. 16|) in the form 

n-1 

5f(T, v(r), ei(r), . . . , £„_i(r)) = -s(r) 5T - p (r) 5v(r) + £ ^ <%(r) (111.17) 

»=i 

where 

s(r) = f(T, v(r), ^(r), . . . , &,_i(r)) = / (T, w(r), &(r), . . . , £„_i(r)) (111.18) 

Eq. (|III.17|) has the form of the usual Gibbs relation for a homogeneous mixture. For an inhomogeneous mixture, 
however, the validity of such a relation is not obvious. Eq. piI.17[) implies that with respect to the variations of the 
thermodynamic variables the specific Helmholtz energy is still homogeneous of the first order. We will call Eq. 1 7|) 
the ordinary Gibbs relation. With the help of Eq. (IIII.10D and Eq. pil,17[) we obtain the Gibbs-Duhem relation 



n-1 

s(r) ST - v(r) 8p (r) + Sfi m + ^ &( r ) = 

i=i 



(111.19) 
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Using the following conditions, which are true for equilibrium 

(111.20) 



VT(r) =0 
VMn(r) =0 



V^i(r) = 0, for i = l,n-l 
V Q cr Q/3 (r) =0 

and Eq. pil.lOj) together with Eq. (fTTTT8]) we obtain 

n-l 

Vp /(r) = -p (r) «(r) + £ ^ &(r) + v(y) V Q 7a/3 (r) (111.21) 

i=i 

We will call Eq. (|III.21[) the spatial Gibbs relation. As the temperature is independent of the position, the expected 
— s(r) V pT term is zero. 

C. Equilibrium surface. 

Away from the surface c(r) — > c, £fc(r) — * £fc an d Vc(r) — > 0, V^(r) — > and we have the homogeneous mixture. 
Thus, the usual thermodynamic relations for the homogeneous mixture are valid. The specific Helmholtz energy is 
given by 

n-l 

/o(T, c, = Mn, o(T, c, + E ^ o( T ' c ' ^ fc - ^ c ' u ( IIL22 ) 

fc=i 

with the specific entropy, pressure and chemical potentials given by 

sa(T, c, £) = /o(T, c, 

Po(T,c,0 =-|j/o(T, c,e) (IIL23) 
^ fc ,o(T, c, =4;fo(T, c, 0, fc= V^T 
such that the variation of the specific Helmholtz energy 

n-l 

df (T, c, = -« (T, c, dT - Po(T, c, dw + £ ^, o(T, c, d£ k (111.24) 

fc=i 

where ipk,o(T, c, £) = fJ.k,o(T, c, |) - / x n , (T, c, £). 

Quantities derived in Subsec. [Ill Aj are converged in a homogeneous limit in a following way 





-> ipk,o{T, c, £) 


= g|-/o(T, c, 




Hn,o{T, C, f) 


n-l 

= /o(r, c, f) - 2 ipi,o(T, c, £)& +p (T, c, £) 

i=l 


p(r) - 


-» Po^, c, f) 


- c 2 |:/o(r, C, 




"> O~a/3,0(T, C, £) 


= Po(r, C, £)<5 a/3 



(111.25) 



We use these limits to determine the meaning of the derived quantities in the interfacial region, where gradients are 
not negligible. 

In equilibrium ipk and /i„ are everywhere constant. Away from the surface they represent the homogeneous chemical 
potentials, which in equilibrium should be constant everywhere, particulary throughout the whole interfacial region. 
Thus it is natural to identify ipk and [i n with the chemical potentials also within the interfacial region. 

Before determine the meaning of p(r) and a a p(r) we have to resolve an ambiguity in the definition of a a p{v). It 
follows from Eq. (|III.7[) that a constant tensor can be added to <T a p(r) without affecting the validity of this equation. 
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In a homogeneous limit this tensor does not vanish, so it should be present in the homogeneous tensor <J a p,o( r ), 
which is proportional to the homogeneous pressure po{T, c, £). The homogeneous pressure po{T, c, £), however, is 
determined unambiguously by the specified equations of state. It follows then that this constant tensor have to be 
equal to zero and cr a p{r) is given by Eq. (|III.8[) unambiguously. 

Using that if>k and \i n are the chemical potentials and Eq. (|III.10|) it is then also natural to identify p(r) given by 
Eq. pil.6j) with a pressure everywhere. This pressure is not constant, however. The tensor cr Q( g(r) can be identified with 
the tensorial pressure. It is known that at the surface one can speak about the "parallel" and the "perpendicular" 
pressure so the pressure reveals tensorial behavior. For a flat surface, when all the properties change in one 
direction, say x, one can show that a xx (r) is the "perpendicular" pressure and p(r) = (T yy (r) = a zz (r) is the "parallel" 
pressure. For curved surfaces such an identification, however, can in general not be made. 

One can also conclude, that the quantity, determined by Eq. (|HI.18|) is the specific entropy of the mixture in the 
interfacial region. It does not have gradient contributions. This is due to the assumption that the coefficients of the 
square gradient contributions are independent of temperature. We refer to van der Waals [H, 0, H| who discussed this 
property. 

We shall also define other thermodynamic potentials in the square gradient model for the interfacial region. Con- 
sidering Eq. pil.lOp and conforming to Eq. (IIII.25[) we define interfacial molar internal energy, enthalpy and Gibbs 
energy densities as follows 

u(r) = /(r) + s(r)T 

h(r) =u(t)+p(t)v(t) (111.26) 
g(r) = /(r) +p(r)v(r) 

It is important to realize that these thermodynamic relations are true in the interfacial region only by definition. One 
can also find support for these definitions in [7[ where they were considered for a simplified one-component system. 

Using Eq. (|III.17|) and Eq. (|III.26j) . for the internal energy at each point in space we then get the ordinary Gibbs 
relation 

n-i 

5u{s(r),v(r),^(r), . . . , f»_i(r)) = TSs(r) - p (r) Sv(r) + £ fa <5&(r) (111.27) 

i=l 

From Eq. (|III.20|) and Eq. (|III.26[) we get the spatial Gibbs relation 

Vp t*(r) = TVp s(r) - p (r) V v(r) + £ Vi V > &(r) + v{r) V a j a0 (r) (111.28) 

i=l 

One can easily write the Gibbs relations for other thermodynamic potentials. 



IV. GRADIENT MODEL FOR THE SPECIFIC VARIABLES PER UNIT OF VOLUME. 

We write the extended Helmholtz energy per unit of volume as 

f v (r) = f$(T,c)+lC v (c, Vc) (IV.l) 

where 

n 

IC v (c, Vc) = - <j Vci(r)-V Cj (r) (IV.2) 

i J = l 

We use c as short notation instead of whole set {ci, . . . , c„} and Vc as short notation instead of {Vci, . . . , Vc„} in 
Sec. [IV] , It should not be confused with the total molar concentration, used in Sec. |III| . 

The coefficients k& are different from those we used in Sec. [Ill] . One can derive the following relations between 
them. 

n—l n—1 

i,j=l i 
n-l 

«i =cE^kV +ckV (IV.3) 
Kij = c 2 k^ 
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where 



— K in C (IV. 4) 



K — K nn 



and 



n 1 n 1 71 

K?- = K + - {Kj + Kj) + — ( Kjj - 2^ K^Q j 3^(^j+Ki<)Q+7 E K4fc c 4 c 4 ( IV - 5 ) 

n 

where c = Cj, 

We will not repeat in details the procedure, given in Sec. |IIIj and will give only the results here. Using the Lagrange 
method we obtain the expressions for the constant Lagrange multipliers, which are equal to the chemical potentials 

ft n 

»k = j-ifo + K?) ~ E V '«* Vc *) ( IV - 6 ) 



dc k 



and the expression for a pressure p: 



n o n 

»=1 1 i,j=l 

which was defined by 

n 

r(r)=E^c,(r)-p(r) (IV.8) 

i=l 

One can derive the same symmetric tensorial pressure a a p(r) as in Eq. pil.8[) . which obeys the relation Eq. (|III.7[) . 
where the symmetric tension tensor 7 Q( 3(r) is given by 

Varying the total Helmholtz energy <5F[T, c(r), Vc(r)] with respect to the variation of the variables we obtain the 
total thermodynamic differential of the specific Helmholtz energy as 

Sf v (T, c, Vc) = ^6T + J2 Mfe Sc k + W-S&" (IV. 10) 

k=l 

where 

n 

S@ v (c, Vc) = E (KijVci)5cj (IV.ll) 

The total Helmholtz energy variation becomes then 

<5F[T, c, Vc] = ^ dr 5T + E ( IV - 12 ) 

since the boundary integral J s dSn s -5& v disappears. Thus, we will interpret the expression in parenthesis as the 
total thermodynamic differential of the specific Helmholtz energy: 

Sf v (T, c, Vc) = ^ ST + E Mfc (IV. 13) 

fc=i 
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We write Eq. (|IV13[) in a form 

n 

5f v (T, Cl (r), ...,c n (r)) = -s v (r)ST + Y,»iSci{r) (IV.14) 

i=l 

which we will call Eq. (jIV. 14|) the ordinary Gibbs relation. Here 

s v (t) = r (T, ci(r), . . . , c„(r)) = fg(T, d(r), . . . , c„(r)) (IV.15) 

is the specific entropy of the mixture. With the help of Eq. (|IV.8j) and Eq. (jiV.14j) we obtain the Gibbs-Duhem 
relation 

n 

a v (T)ST-Sp(T) + ^2c k (T)Sii k = (IV. 16) 

fe=i 

Using the following conditions, which are true for equilibrium 

VT(r) =0 

V fi k (r) = 0, for k — T~n (IV.17) 
V Q CT Q /3(r) = 

and Eq. pV.8|) together with Eq. 8[) we obtain the spatial Gibbs relation 

n 

V/j f(r) = W V/j Ci(r) + V Q 7o/J (r) (IV.18) 
»=i 

The interfacial specific internal energy, enthalpy and Gibbs energy densities are 

u v {r) = f v {r) + s v (r)T 

n 

h v (r) = s v {r)T+J2f H c i {r) (IV.19) 

i=l 

n 

9 V (*) = E fMCi(r) 

i=l 

The Gibbs relations for the internal energy are 

n 

Su v (s(r), Cl (r), ...,c n (r))=T5 S v (r)+Y / ^Sci(r) (IV.20) 

i=l 

and 

V Q 7 a /3(r) (IV. 21) 

V. TYPICAL PROFILES FOR THE BINARY MIXTURE. 

In order to illustrate the results, which one can obtain using the above procedure, we have applied it to a special 
case. This requires a number of approximations, connected with the specific mixture and the geometry. We consider 
a flat liquid-vapor interface of the binary mixture of cyclohexane (1 st component) and n-hexane (2 nd component). 
The equilibrium profiles are obtained using Eq. (|III.5[) . The values for the coexisting chemical potentials are obtained 
from the coexisting conditions using the van-der-Waals equation of state. For the binary mixture this requires two 
independent thermodynamic properties, which we have chosen to be the equilibrium temperature T = 330 K and 
molar fraction of the first component in the liquid phase = 0.5. As cyclohexane and n-hexane are rather similar, the 
gradient coefficients k\i, k% 2 an d K u were chosen to be equal. We use the value 12E-18 J*m 5 /mol 2 which reproduces 
the experimental value of the surface tension 0.027 N/m. The resulting profiles for the molar concentration, the mole 
fraction of the first component and the tension tensor component "f xx , the integral of which gives the surface tension, 
are given in Figs. [TJS1 This paper is the first in a sequence of papers. In the second paper we will extend the present 
analysis to non-equilibrium mixtures. Details of the numerical procedure will be given in a third paper which discusses 
the validity of local equilibrium for the Gibbs surface. 




FIG. 1: Molar concentration profile 
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FIG. 2: Molar fraction profile 
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FIG. 3: Tension j xx profile 

VI. DISCUSSION AND CONCLUSIONS. 

In this paper we have established the framework of the gradient model for the liquid-vapor (or, alternatively, 
liquid-liquid) interface in an isotropic non-polarizable mixture. It is necessary that the homogeneous Hclmholtz 
energy /q allows solutions which imply equilibrium coexistence between different phases. Otherwise we only have the 
homogeneous phase. Standard mixture theories [12| give a Helmholtz energy which allows liquid-vapor coexistence. 

Using the assumption that in the interfacial region the fluid can be described by the local densities and their 
gradients, we have extended the gradient models, used to describe one-component fluids and binary mixtures, to three- 
dimensional multi-component mixtures. The condition which the system should satisfy is that the total Hclmholtz 
energy is minimal. With the help of Lagrange method is was possible to derive the equations, which the profile 
distribution should satisfy, given the fixed total molar content of the components. The Lagrange multipliers are equal 
to the chemical potentials of the coexisting liquid and vapor. It was also possible to determine the pressure behavior 
in the interfacial region. It is crucial, that the pressure has a tensorial behavior. The difference between the tensorial 
part of the pressure tensor and the scalar part determines the surface tension. 

An important observation is the ambiguity in the determination of the local thermodynamic potentials, for instance 
the specific Helmholtz energy. While the total Hclmholtz energy is unique and has a minimum, the specific Hclmholtz 
energy is not unique. One can add a term which is the divergence of some vector field, if the normal component of 
this field vanishes on the system boundaries, without affecting the total Helmholtz energy. This general observation 
in the context of the gradient model implies that the density gradients are taken equal to zero on the boundary. It 
must be emphasized, that for realistic boundaries these gradients are not zero. We refer to the wall-theorem [1 31 ] in 
this context. We take the freedom to set them equal to zero assuming, that the boundary layer does not affect the 
properties of the interface we want to study. In the Appendix [X] we show that this ambiguity does not affect the 
results. 

As one can see from Eq. 10|) and Eq. (|IV.8[) . it is possible to relate the thermodynamic variables for an 
inhomogeneous fluid in the same way as it is done for a homogeneous one. However, unlike the homogeneous mixture, 
these variables contain gradient contributions. This means that local equilibrium for such a system is not satisfied. 
The local behavior of the mixture is determined not only by it's local properties but also by it's nearest surroundings. 
Moreover, in contrast to the homogeneous description, the local properties now vary in the space. 
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We have given explicit expressions for each thermodynamic quantity in the interfacial region. We have also deter- 
mined how the thermodynamic potentials change with the change of the variables they depend on. An important 
part of the thermodynamic description is the relations between the rate of change of the thermodynamic variables, 
the Gibbs relations. In contrast to a homogeneous system, for an inhomogeneous system, in particular the interfacial 
region, thermodynamic variables vary also in space. Thus one can speak about the relation between the rates of 
change of thermodynamic variables for a given point in space, the ordinary Gibbs relation. One has also to speak 
about the rates of change of the thermodynamic variables in space, the spacial Gibbs relation. Even though the 
thermodynamic potentials, particularly the specific Helmholtz energy, depend on the spatial derivatives of the den- 
sities, we have shown that variation of these gradients do not contribute to the ordinary Gibbs relations. Thus, the 
ordinary Gibbs relations have the ordinary form of the Gibbs relations for the homogeneous mixture. The important 
observation here is, however, that the ordinary Gibbs relations relate the inhomogeneous thermodynamic variables, 
i.e. those, which contain the gradient contribution. As the spatial derivatives of the temperature, chemical potentials 
and pressure tensor are zero in equilibrium, we can determine the spatial Gibbs relation. The new term which appears 
because of the inhomogcncity is V Q 7^/3 (r), which is only unequal to zero close to the surface. 

For temperatures far from the critical point the surface thickness is known to be very small (in the sub-nm range) . 
This imposes an upper limit to the values of the coefficients k, /tj, Kij and 

In previous sections we used different specific variables. For molar densities we had the coefficients k, Kj and 
and for densities per unit of volume k^. One can determine the relations between these coefficients and verify, that 
all the quantities, determined in Sec. |IIIj and Sec. |IV| are the same. Thus, Eq. (|III.6[) and Eq. (|IV.7|) give the same 
quantity p (r), Eq. (|III.9[) and Eq. pV.9|) - the same 7 Qj g(r). And <r aj g(r), which is given by Eq. (|III.8[) is the same for 
both sets of variables. \x n in Eq. (|III.5|) and Eq. pV.6p is the same and ipk taken from Eq. (|III.5|) are equal to /Uj. — fi n 
taken from Eq. (|IV.6[) . This shows that the inhomogeneous equilibrium description is independent of the choice of 
independent variables. This is similar to the description of the homogeneous equilibrium phase. 

The analysis in this paper gives the basis to extend the description to non-equilibrium systems. For one-component 
systems, in which the properties varied only in one direction, such an extension was given by Bedeaux et. al. 
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APPENDIX A: ON THE AMBIGUITY IN THE SPECIFIC QUANTITIES. 

We show here that ambiguities present in the definition of the specific Helmholtz energy and the total thermo- 
dynamic differential of the specific Helmholtz energy do not affect the validity of all the thermodynamic relations, 
derived in the present article. We do this for molar specific variables, the same arguments can be used for the specific 
variables per unit of volume. 

The total Helmholtz energy of a mixture J y dr c(r)/(r) can be expanded around homogeneous state as 



F = f v dr c/ (T, c, + 



n — 1 . , n—1 , . 



2 ^ — * ••1] 



(A.l) 



-k^ v 2 c + e 4 2 ] v 2 e 



where the series is truncated after the second order terms. As it is shown in Eq. (|II.2[) this can be rearranged as 
following: 



F = J v dr cf (T,c,0- 



n—1 n—1 

+±k|Vc| 2 + E KiVc-VCi + l E ««V&-V& 

i—1 i,j=l 
n—1 

+v-(^ 2 )vc+ e «rvfi) + .. 
i=i 



(A.2) 
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where 



K 



(2) 



(i) _ a K < 2) _ £« 



(2) 



are the coefficients used in Eq. (|III.2I) . 

Let us define 

n— 1 7i — 1 

/i(r) = c/ (T, c, £) + |/s |Vc | 2 + E «i Vc-VC + |E «tf 

Z— 1 2,J=1 
ft— 1 fo \ 

i=l 

so that F = f v drc(r)fi(r) = J v dr c(r)/ 2 (r), since J 7 drc(r)(/ 2 (r) — /i(r)) = 0. One can also define 

n— 1 / n— 1 

pi(r) =/i„+ E -/i(r) =c 2 |;(/o + /C) -cV-( K Vc+ E «i V£ 

<=1 v i=l 

ri — 1 n— 1 . . 

p 2 (r) = /i„ + E V^W - / 2 (r) =pi(r) - V-( K ( 2 >Vc + E «r V &) 

i=l j=l 



(A.4) 



(A.5) 



and follow the same procedure as in Subsec. [Ill Aj . (Note, that the second order terms should be added to the 
Euler-Lagrange equations in thus case). Then, as it follows from the known theorem of the variational calculus, one 
obtains for the chemical potentials /z„ and 



2—1 i=l 

9 1 71—1 

^fe = ^ (/i) - - V- (k* Vc + K lfc V&J , k = T~^l 



(A.6) 



and for the pressure tensor a a fj{v) 

Ca/3(r) = pi(r) S a /3 + 7a/3(r) (A. 7) 

for the both definitions Pi(r) and P2( r )- As one can see, there is no ambiguity in the expressions for the chemical 
potentials and the pressure tensor. They are determined irrespectively of the ambiguity in the definition of the specific 
Helmholtz energy. Thus, it is natural to use /i(r) as the specific Helmholtz energy and use pi(r) as a pressure. As we 
could see in Subsec. C| only p±(r) has a physical meaning but not P2( r )- With such a choice all the thermodynamic 
quantities, derived in the paper do not contain any ambiguity. 



APPENDIX B: SYMBOLS LIST AND NOTATION CONVENTION. 



- contraction sign 

l,n enumeration of all integers from 1 to n 

a, f3 - cartesian indices 

8 total thermodynamic differential 

S a /3 Kroneker symbol 



V [ 1/m ] - nabla operator 

V Q [ 1/m ] - partial derivative with respect to x a 

7 a ^ [ Pa ] - tension tensor 

A [ 1/m 2 ] - Laplace operator 

K, [ J/mol ] - square gradient contribution 

fC v [ J/m 3 ] -»--»- 
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K 


[ J m 5 /mol 2 ] 


- 


square gradient coefficients for molar units 


Ki 


[ J m/mol ] 
[ J/m 3 ] 


- 


>> >> 


Kij 


- 


>> >> 


V 

ij 


[ J m 5 /mol 2 ] 


- 


square gradient coefficients for volume units 


Hk 


[ J/mol ] 


- 


molar chemical potential of component k 


Hk 


[J/kg], [J/mol] 


- 


matter chemical potential of component k 


Pk,0 


[ J/mol ] 


- 


homogeneous molar chemical potential of component k 




[ J/kg ], [ J/mol ] 


- 


homogeneous matter chemical potential of component k 


V 


[ mol ] 


- 


total number of moles 


l'k 


[ mol ] 


- 


number of moles of component k 


p 


[ kg/m 3 ] 


- 


total mass density 


Pk 


[ kg/m 3 ] 


- 


mass density of component k 




[ Pa ] 


- 


pressure tensor 




["] 


- 


set of molar fractions {£1, . . . , £ n -i} 




["] 


- 


molar fraction of component k 


& 


["] 


- 


matter fraction of component k 




[ J/mol ] 


- 


molar reduced chemical potential of component k 




[J/kg], [J/mol] 


- 


matter reduced chemical potential of component k 


1pk,0 


[ J/mol ] 


- 


homogeneous molar reduced chemical potential of component k 




[ J/kg ], [ J/mol ] 


- 


homogeneous matter reduced chemical potential of component k 


c 


[ mol/m 3 ] 


- 


total molar concentration 


c 


[ mol/m 3 j 


- 


short notation for the set {ci, . . . , c n } of molar concentrations (only in Sec. 


Ck 


[ mol/m 3 ] 


- 


molar concentration of component k 


F 


[J] 


- 


total Helmholtz energy 


f 


[ J/mol ] 


- 


specific Helmholtz energy per mole 


fo 


[ J/mol ] 


- 


homogeneous specific Helmholtz energy per mole 


r 


[ J/m 3 ] 


- 


specific Helmholtz energy per unit of volume 


Jo 


[J/kg], [J/mol] 


- 


homogeneous specific Helmholtz energy per unit of matter 


Jo 


[ J/m 3 ] 


- 


homogeneous specific Helmholtz energy per unit of volume 


i,j,k 


["] 


- 


component number 


m 


[kg] 


- 


total mass 


mk 


[kg] 


- 


total mass of component k 


m 


[ kg ], [ mol ] 


- 


total matter amount 


Ul k 


[ kg ], [ mol ] 


- 


total matter amount of component k 


ml 


[ kg/m 3 ], [ mol/m 3 ] 


- 


matter density of component k: amount of component k per unit of volume 


n 


["] 


- 


number of components 


P 


[ Pa ] 


- 


scalar pressure 


Po 


[ Pa ] 


- 


homogeneous pressure 


r 


[m] 


- 


position 


S 


[J/K] 

[m 2 ] 


- 


total entropy 


S 


- 


the boundary surface 


s 


[ J/(K mol) ] 


- 


specific entropy per mole 


s v 


[ J/(K m 3 ) ] 


- 


specific entropy per unit of volume 




[J/(K kg)], [J/(K mol)] 


- 


homogeneous specific entropy per unit of matter 


So 


[ J/(K m 3 ) ] 


- 


homogeneous specific entropy per unit of volume 


T 


[K] 


- 


temperature 


U 


[J] 


- 


total internal energy 


u 


[ J/mol ] 
[ J/m 3 ] 


- 


specific internal energy per mole 


V 

a 


- 


specific internal energy per unit of volume 


u? 


[J/kg], [J/mol] 
[ J/m 3 ] 


- 


homogeneous specific internal energy per unit of matter 


ul 


- 


homogeneous specific internal energy per unit of volume 


V 


[ m d ] 




total volume 


V 


[ m 3 /mol ] 




specific volume per mole 




[ m 3 /kg ] 




specific volume per unit of mass 


v m 


[ m 3 /kg ], [ m 3 /mol ] 




specific volume per unit of matter 




[ m ] 




cartesian coordinate in direction a 
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